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Abstract 

Resource  envelopes  provide  the  tightest  exact  bounds  on  the 
resource  consumption  and  production  caused  by  all  possible 
executions  of  a temporally  flexible  plan.  We  present  a new 
class  of  algorithms  that  computes  an  envelope  in 
0(Maxflow(n,  m,  U))  where  n,  m and  U measure  the  size  of 
the  flexible  plan.  Tliis  is  an  0(n)  improvement  on  the  best 
complexity  bound  for  an  envelope  algorithm  known  so  far 
and  makes  envelopes  more  amenable  to  practical  use  in 
scheduling.  The  reduction  in  complexity  depends  on  the  fact 
that  when  the  algorithm  computes  the  constant  segment  i of 
the  envelope  it  makes  full  reuse  of  the  maximum  flow  used 
to  obtain  segment  /-/. 


Resource  Envelopes 

The  execution  of  plans  greatly  benefits  from  temporal 
flexibility.  Fixed-time  plans  are  brittle  and  may  require 
extensive  replanning  due  to  execution  uncertainty. 
Moreover,  when  plans  must  deal  with  uncontrollable 
exogenous  events  (Morris  et  al.,  2001)  temporal  flexibility 
cannot  be  avoided.  However,  effective  algorithms  to  build 
temporally  flexible  plans  are  rare,  especially  when 
activities  produce  or  consume  variable  amounts  of  resource 
capacity.  A major  obstacle  is  the  difficulty  of  assessing  the 
resource  needs  across  all  possible  plan  executions. 

Methods  are  available  to  compute  resource  consumption 
bounds  (Laborie,  2001;  Muscettola,  2002).  In  particular, 
(Muscettola,  2002)  proposes  a polynomial  algorithm  to 
compute  a resource  envelope,  the  tightest  of  these  bounds. 
By  being  the  tightest,  resource  envelopes  can  potentially 
save  an  exponential  amount  of  search  (through  early 
backtracking  and  solution  detection)  when  compared  to 
using  looser  bounds.  Also,  methods  that  compute  resource 
envelopes  identify  maximally  matched  sets  of  resource 
consumer/producers  that  balance  each  other  for  any  plan 
execution.  This  and  other  structural  information  could  be 
crucial  in  minimizing  the  search  space  and  suggesting 
effective  scheduling  heuristics,  potentially  enabling  new 
classes  of  highly  efficient  schedulers. 
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However,  preliminary  comparative  studies  of  scheduling 
algorithms  using  envelopes  appear  not  to  show  a 
computational  advantage  with  respect  to  using  more 
traditional  heuristic  methods  based  on  fixed-time  resource 
profiles  (Pollicella  et  al.,  2003).  Since  computing  envelopes 
is  more  computational  expensive  than  building  a fixed-time 
profile,  it  is  critical  to  ensure  that  the  balance  between 
computation  cost  and  increased  structural  information 
extracted  form  the  envelope  is  advantageous.  Making  the 
trade-off  advantageous  requires  two  complementary 
approaches.  The  first  reduces  the  cost  of  computing  an 
envelope;  the  second  devises  new  envelope  analysis 
methods  to  extract  useful  heuristics. 

In  this  paper  we  address  the  problem  of  cost  reduction. 
Currently,  the  resource  envelope  algorithm  known  to  have 
the  best  asymptotic  complexity  (Muscettola,  2002) 
computes  all  piecewise-constant  segments  of  the  envelope 
through  as  many  as  2n  stages,  where  n is  the  number  of 
events  (start  or  end  of  activities)  in  the  flexible  plan.  Each 
stage  computes  a maximum  flow  and  therefore  the  overall 
complexity  of  the  method  is  0(n  Maxflow(n,m,U))  where  m 
is  the  number  of  temporal  constraints  between  activities  in 
the  plan,  U is  the  maximum  level  of  resource  production  or 
consumption  at  some  activity,  and  Maxflow(n,  m,  f/j  is  the 
asymptotic  cost  of  the  maximum  flow  algorithm. 

This  staged  method,  however,  can  be  significantly 
improved  since  at  each  stage  a full  maximum  flow  for  the 
entire  flexible  plan  is  recomputed  from  scratch.  Cost 
reduction  could  be  obtained  through  an  incremental  flow 
method.  Starting  from  the  maximum  flow  at  one  stage,  the 
maximum  flow  for  the  next  stage  is  obtained  by  minimally 
reducing  flow  when  deleting  nodes  and  edges,  and  by 
minimally  increasing  flow  when  adding  new  nodes  and 
edges  (Kumar,  2003).  However,  without  appropriately 
ordering  flow  reductions  and  increases,  the  asymptotic 
complexity  may  not  improve  (at  it  appears  to  be  the  case  in 
(Kumar,  2003)). 

In  this  paper  we  introduce  an  incremental  method  that 
provably  computes  an  envelope  in  0(Maxflow(n,  m,  U))  for 
a large  class  of  maximum  flow  algorithms.  This  reduction 
of  complexity  is  significant.  Experimental  analysis  has 
shown  that  the  practical  cost  of  maximum  flow  is  usually 
as  low  as  0(n  (Ahuja  et  al,  1993).  This  compares  well 
with  0(n  log  n),  the  cost  of  building  resource  profiles  for 
fixed  time  schedules. 
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This  paper  is  organized  as  follows.  We  first  give  a 
succinct  introduction  to  the  resource  envelope  problem  and 
the  staged  envelope  algorithm  in  (Muscettola,  2002).  Next 
we  present  the  new  incremental  algorithm  and  identify  all 
sources  of  performance  improvements.  We  then  prove  the 
complexity  result,  discuss  implementation  improvements 
when  using  preflow-push  algorithms  and  conclude  by 
discussing  future  work. 

Staged  Computation  of  Envelopes 

In  this  section  we  outline  the  envelope  problem  and  the 
staged  algorithm  that  solves  it.  For  a complete  discussion, 
see  (Muscettola,  2002). 

Figure  1 shows  an  activity  network  with  resource 
allocations.  The  network  has  two  time  variables  per 
activity,  a start  event  and  an  end  event  (e.g.,  and  for 
activity  A^),  a non-negative  flexible  activity  duration  link 
(e.g.,  [2,  5]  for  activity  A),  and  flexible  separation  links 
between  events  (e.g.,  [0,  4]  from  to  Two  additional 
events  and  define  a time  horizon  within  which  all 
events  occur. 

Time  origin,  events  and  links  constitute  a Simple 
Temporal  Network.  To  describe  resource  production  and 
consumption  each  event  also  has  an  allocation  value  r(e) 
(e.g.,  r(ej  = -2),  a numeric  weight  that  represents  the 
amount  of  resource  allocated  when  the  event  occurs.  We 
will  assume  that  all  allocations  refer  to  a single,  multi- 
capacity resource.  The  extension  to  multiple  resources  is 
straightforward.  If  the  allocation  is  negative  an  event  e~  is  a 
consumer,  if  it  is  positive  e^  is  a producer.  We  assume  that 
the  temporal  constraints  are  consistent  which  means  that 
for  any  pair  of  events  the  shortest  path  \eje^\  from  e^  to  e^  is 
well  defined.  Each  event  e can  occur  within  its  time  bound, 
between  the  earliest  time  et(e)  = and  the  latest  time 
lt(e)  = \Tje\,  The  triangular  inequality  \e^e^\  + \e^e^\ 

holds  for  any  three  events  e^  and  e^ 

<e2s,  3>  [2,3]  <C2c,  2> 


Figure  1:  An  activity  network  with  resource  allocations 

Informally,  a flexible  plan  is  resource  consistent  if  the 
duration  and  separation  links  induce  appropriate  necessary 
precedence  relations  between  consumers  and  producers. 
These  relations  should  guarantee  that,  when  a consumer 
occurs,  the  total  resource  level  due  to  consumers  and 
producers  that  cannot  occur  after  it  must  be  at  least  as  high 


as  the  new  consumption.  A similar  condition  applies  to  the 
correct  occurrence  of  a producer.  The  full  information  on 
the  necessary  precedence  relations  is  captured  the  anti- 
precedence graph  Aprec,  a graph  that  contains  a path 
between  any  two  events  Cj  and  e^  if  and  only  if  \ej  cj  < 0. 
Figure  2 depicts  an  anti-precedence  graph  of  the  network  in 
Figure  1 with  each  event  labeled  with  its  time  bound  and 
resource  allocation.  We  use  anti-precedence  graphs  rather 
than  the  most  customary  precedence  graphs  (Laborie,  2001) 
to  simplify  the  construction  of  the  auxiliary  maximum  flow 
problem  that,  as  we  will  see,  is  fundamental  for  the 
computation  of  envelopes. 

We  can  now  formally  define  a resource  envelope.  For 
any  subset  of  events  A,  the  resource  level  increment  of  A is 
A(A)  = 0 if  A = 0.  and  A(A)  = if  A If  5 is  the 

set  of  all  possible  consistent  time  instantiations  for  all 
events  and  r is  a time  within  the  time  horizon,  the  resource 
level  at  time  t for  a specific  time  instantiation  .y  e 5 is  L/t) 
= A(Ej(t)),  Here  Ej(t)  is  the  set  of  events  e which  occur  at 
or  before  t in  s.  The  maximum  resource  envelope  is  L^Jt) 
= max^^^L/t)  and  the  minimum  resource  envelope  is  L^Jt) 


<[4,  10],  3>  <[6,  13],  2> 


Figure  2:  Anti -precedence  graph  with  time  bounds  and 
resource  allocations 

= min^^gL/t),  Since  can  be  computed  with  obvious  term 
substitution  on  the  method  that  computes  we  only 
focus  on  L^, 

To  compute  the  resource  envelope  at  time  t we  partition 
all  events  into  three  sets  depending  on  the  position  of  their 
time  bound  relative  to  r:  1)  the  closed  events  C,  that  must 
occur  before  or  at  r,  i.e.,  such  that  that  lt(e)  < t\  2)  the 
pending  events  /?,  that  can  occur  before,  at  or  after  r,  i.e., 
such  that  et(e)  < t < lt(e)\  and  3)  the  open  events  O,  that 
must  occur  strictly  after  r,  i.e.,  such  that  et(e)  > t. 

Any  resource  level  increment  L/r)  will  always  include 
the  contribution  of  all  events  in  C,  and  none  of  those  in 
but  may  include  only  some  subset  of  events  in  i.e.,  only 
those  that  are  scheduled  before  t in  s.  It  is  possible  to  show 
that  this  subset  must  be  a predecessor  set  such  that  if 
eeP  and  e*  follows  e in  Aprecy  then  e'sP,  We  call  PmcJPt) 
the  (possibly  empty)  predecessor  set  with  maximum  non- 
negative resource  level  increment. 

The  fundamental  result  reported  in  (Muscettola,  2002)  is 
that  LfruJt)  can  be  determined  from  the  following  equation. 

Equation  1:  L^(t)  = A(Q^A(P„JR)) 
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Figure  3:  A resource  increment  flow  network 

Assuming  that  time  bound  information  is  available  for 
all  events,  the  computation  cost  of  A{C)  is  0(n),  The  cost 
of  computing  P„JR)  determines  the  asymptotic  cost  of 
L,nax(0-  We  can  compute  P^JR)  by  solving  a maximum 
flow  problem  on  an  auxiliary  flow  network  F (R),  the 
resource  increment  flow  network  for  R^. 

The  formal  definition  of  a resource  increment  flow 
network  can  be  found  in  (Muscettola,  2002).  As  an 
example.  Figure  2 gives  F(RJ  for  the  activity  network  in 
Figure  1.  The  network  has  a node  for  each  event  in  R^,  an 
infinite  capacity  flow  edge  between  two  events  for  each 
edge  in  Aprec  (see  Figure  2),  an  edge  from  the  source  crto 
a producer  with  capacity  equal  to  the  producer’s  allocation, 
and  an  edge  from  a consumer  to  the  sink  r with  capacity 
equal  to  the  opposite  of  the  consumer’s  allocation. 

A complete  discussion  of  maximum  flow  algorithms  can 
be  found  in  (Cormen,  Leiserson  and  Rivest,  1990).  Here  we 
only  highlight  a few  concepts  that  we  will  use  in  the 
following.  A flow  is  a function  e^)  of  pair  of  events  in 
F{R)  that  is  skew-symmetric,  \,t,,f[e^,  ej  = ef),  for 
each  edge  e^—>e^  has  a value  no  greater  than  the  edge’s 
capacity  c{e^,  ef)  (assuming  capacity  zero  if  the  edge  is  not 
in  F(Rf),  and  is  balanced,  i.e.,  the  sum  of  all  flows  entering 
an  event  must  be  zero.  A pre-flow  is  a function  defined 
similarly  but  that  relaxes  the  balance  constraint  by  allowing 
the  sum  of  pre-flows  entering  a node  to  be  positive.  The 
total  network  flow  is  defined  as  e)  = 

The  maximum  flow  of  a network  is  a flow  function  f^ 
such  that  the  total  network  flow  is  maximum. 

A fundamental  concept  in  the  theory  of  flows  is  the 
residual  network  for  a particular  flow,  a graph  with  an 
edge  for  each  pair  of  nodes  in  F(R)  with  positive  residual 
capacity,  i.e.,  the  difference  c(Cj,  ej  - fie j,  efj  between 
edge  capacity  and  flow.  Each  residual  network  edge  has 
capacity  equal  to  the  residual  capacity.  An  augmenting  path 
is  a path  connecting  o*  to  r in  the  residual  network.  The 
existence  of  an  augmenting  path  indicates  that  additional 
flow  can  be  pushed  from  a to  r.  Alternatively,  the  lack  of 
an  augmenting  path  indicates  that  a flow  is  maximum. 

We  can  compute  P„JR)  according  to  the  next  theorem. 

Theorem  2:  [Theorem  1 in  (Muscettola  2(X)2)]  P^(R) 
is  the  (possibly  empty)  set  of  events  that  are  reachable  from 
the  source  a in  the  residual  network  of  some  ofF(R). 


From  Equation  1 and  Theorem  2 (Muscettola,  2002) 
derives  a staged  envelope  algorithm  as  follows.  Consider  a 
time  corresponding  either  to  the  earliest  or  the  latest  time 
of  some  event.  In  a network  with  n events  there  are  at  most 
2n  such  times.  Since  the  envelope  level  can  only  change  at 
one  of  these  times,  the  algorithm  computes  a different  level 
for  each  of  them.  At  a particular  t^  the  algorithm  determines 
the  corresponding  closed  event  set  C.  and  pending  event  set 
/?;,  builds  F(R),  computes  one  of  its  maximum  flow  using 
some  appropriate  maximum  flow  algorithm,  determines 
PmJR)  according  to  Theorem  2,  and  computes  L^(t) 
according  to  Equation  1.  It  is  easy  to  see  that  the  worst-case 
time  complexity  of  this  algorithm  is  0(n  Maxflow(n,  m,  U)) 
where  Maxflow(n,  m,  U)  is  the  worst  time  complexity  of 
the  maximum  flow  algorithm  used. 

Incremental  Computation  of  Envelopes 

In  the  envelope  algorithm  previously  described,  maximum 
flows  are  recomputed  from  scratch  for  each  F(RJ.  Assume 
that  the  times  t.  are  sorted  in  order  of  increasing  value.  To 
reduce  the  cost  of  computing  the  maximum  flow  for  F(Rf, 
we  will  follow  the  approach  of  reusing  as  much  as  possible 
of  the  maximum  flow  computed  for  F (R.J.  Our  theory  is 
independent  from  specific  maximum  flow  algorithms. 
Instead  we  build  our  argument  on  general  properties  of  the 
resource  increment  flow  networks. 

As  in  (Muscettola,  2002)  the  fast  envelope  algorithm 
operates  by  the  modular  application  of  a maximum  flow 
algorithm  to  a well  defined  sequence  of  resource  increment 
flow  network.  The  generality  of  the  following  theory  opens 
the  possibility  of  studying  which  flow  algorithm  could  be 
most  appropriate  for  different  kinds  of  plan  topologies. 

Sources  of  Incremental  Envelope  Speedup 


Figure  4;  Incremental  modification  of  a resource  flow  network 

At  time  t.  the  set  of  pending  events  can  undergo  two 
modifications.  First,  the  events  SC.  = - R.  move  from 

R.  j to  C^.  These  are  events  e such  that  r.  = lt(e).  Second,  the 
events  = R.  - R^  j move  from  0.  ^ to  R.,  These  are  the 
events  e such  that  t.  = et(e).  For  example,  consider  the 
activity  network  in  Figure  1 and  the  process  through  which 
R.  J for  t^  j-  3 is  transformed  into  R.  for  t.  = 4.  This  is 
described  in  Figure  4 where  at  time  4 the  grayed  part  of  the 
network  is  deleted  and  the  emphasized  part  of  the  network 
is  added.  In  particular,  we  have  SC^  = (ej  and  SR.  = (ej. 
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For  completeness,  we  note  that  F(SC.)  consists  of  node 
and  edge  while  F ( SR J consists  of  node  ^2,  and  edge 
cr-^2s*  other  added  and  deleted  edges  are  connectives 
from  F(R.,- SCJ  to  F( SC)  (edges  e.—^.  and  and 

from  F(SR)  and  F(R^_j-  SC)  (edge  ^2.”^/^)- 

The  sets  SC^  and  satisfy  the  following  fundamental 
properties. 

Lemma  3:  SC^  is  a predecessor  set  contained  in  R,  j.  SR^  is 
the  complement  of  predecessor  set  R.j  in  R^. 

Proof:  We  only  give  the  proof  for  SC.  since  the  one  for 
is  analogous.  We  need  to  show  that  no  predecessor  of  an 
event  in  SC^  can  belong  to  its  complement  in  R.  j,  In  other 
words,  given  e,  e S(C)  and  e^  € R^j-SC^,  using  the 
definition  of  anti-precedence  graph  and  predecessor  set,  it 
must  be  {e^e^l  > 0.  From  the  definition  of  SC.  we  have  lt(e) 
= t.  and  Idej  >t.+J\  From  the  triangular  inequality  applied 
to  the  latest  times  of  Cj  and  ^2!’ 

deduce  \e^e2\  ^IHe^)  -lt(e)  >t.-hJ-t.  = l>  O.D 

Lemma  3 determines  what  flow  edges  are  eliminated 
when  SC.  is  deleted  and  what  are  added  when  is  added. 
In  particular,  we  can  only  delete  edges  that  enter  events  in 
SC^  or  go  from  SC.  to  r.  Similarly,  we  can  only  add  edges 
that  exit  events  in  <5??.  or  go  from  a to  Unlike  previous 
proposals  for  incremental  envelope  calculation  (Kumar, 
2003),  our  method  relies  on  events  and  edges  exiting  and 
entering  the  current  flow  network  in  a well  defined  order. 
This  is  the  primary  key  to  reducing  complexity. 

Directly  related  to  Lemma  3 is  the  possibility  of 
computing  the  maximum  flow  of  F (R)  by  incrementally 
modifying  the  fiow  of  F(Ri.i),  reusing  both  flow  values  and 
intermediate  data  structures  across  successive  invocations 
of  a maximum  flow  algorithm.  We  will  prove  that  our  flow 
modification  operators  guarantee  the  maximality  of  each 
intermediate  flow.  Maintaining  intermediate  flow 
maximality  and  reusing  data  structures  are  keys  to  reduce 
complexity  for  different  kinds  of  maximum  flow 
algorithms. 

A final  factor  is  minimizing  the  size  of  each  intermediate 
flow  network.  We  will  show  that  as  soon  as  the  weight  of 
an  intermediate  is  used  in  the  envelope  calculation, 
F (Pmax)  and  all  of  its  connecting  edges  can  be  safely 
eliminated  from  further  consideration.  This  reduces  flow 
network  size  and  further  contributes  to  cost  reduction. 

The  argument  to  construct  the  fast  envelope  algorithm 
will  proceed  as  follows.  On  the  basis  of  Lemma  3,  we  first 
define  flow  transformation  operators  that  apply  to  network 
additions/deletions  occurring  between  time  /,./  and  time  r,. 
The  operators  define  the  sequence  of  maximum-flow 
problems  that  need  to  be  solved.  Then  we  identify  a flow 
separation  property  that,  once  applied  to  the  sequence  of 
maximum  flow  problems,  further  reduces  the  size  of  each 
step’s  maximum  flow  problem.  Finally,  we  determine  a 


* For  ease  of  exposition  we  assume  discrete  time  although 
the  theory  applies  also  to  continuous  time  with  appropriate 
modifications. 


recursive  equation  that  computes  LmeJU)  as  a function  of 
FmaJ^hi)  and  the  weights  of  event  sets  deduced  from  the 
application  of  flow  transformation  operators. 

Flow  Modification  Networks 

The  philosophy  of  each  flow  transformation  operator  is 
similar  to  that  used  by  the  flow  augmentation  method  in 
maximum  flow  theory.  However,  we  use  this  method  more 
generally  not  only  to  augment  flow  but  also  to  shift  flow 
around  the  network  and  to  reduce  flow.  The  general  idea  is 
the  following.  Given  a flow  network  F and  one  of  its 
maximum  flows  /,  an  operator  first  defines  an  auxiliary 
flow  transformation  network  then  finds  one  of  its 
maximum  flows  and  finally  produces  a flow/,^^  =/+/r. 
Each  Ft  consists  of  selected  edges  in  the  residual  network 
of  Ffor  /.  Since  the  properties  of  flows  are  preserved  in  the 
sum  of  a flow  of  Fand  a flow  of  its  residual  network,  f^^  is 
also  a flow  for  network  F. 

Consider  now  the  resource  increment  flow  network 
F (R.J  at  stage  i-1  and  assume  that  the  set  of  new  closed 
events  SC.  is  not  empty.  At  stage  i all  events  in  SC.  and  all 
of  its  incoming  and  outgoing  edges  will  be  deleted.  This 
also  means  that  any  flow  that  at  the  end  of  stage  i-I  enters 
SC.  will  necessarily  have  to  be  zeroed,  i.e.,  pushed  back 
into  F (RiJ.  The  value  of  this  flow  is  the  sum  of  the 
residual  capacities  of  all  edges  e^-^^  where  ej  e SC.  and 
e2ER.  j - SC..  When  pushed  back,  this  flow  can  follow  two 
routes.  The  first  reaches  r through  some  non-saturated 
exiting  edges  of  F (R. -SC).  If  after  having  followed  the 
first  route  some  flow  is  still  flowing  on  some  Cj—^2 
flow  cannot  reach  x any  more,  a second  route  allows 
reversing  the  remaining  flow  all  the  way  to  <7.  We  call  this 
flow  push-back  operation  eiflow  contraction.  The  first  flow 
route  corresponds  to  a flow  shift  and  the  second  one  to  a 
flow  reduction.  For  example,  consider  the  network  in  Figure 
4.  Assume  that  at  t^3  it  is/„Je,^  t)  = 4,  f^ie^^  t)  = 1 and 
t)  = 2.  At  t=4  the  elimination  of  e^^  requires  pushing 
back  4 units  of  flow.  Three  of  these  units  can  still  reach  r 
by  being  shifted  to  e^^—>T.  Only  one  unit  of  flow  needs  to  be 
pushed  back  to  <7,  If  we  pushed  four  units  of  flow  back  to  <j 
without  shifting  (as  in  (Kumar,  2003)),  later  we  would  need 
to  push  again  three  units  of  flow  from  erto  rto  ensure  flow 
maximality.  This  repetition  of  work  affects  worst-case 
asymptotic  complexity. 

Assume  now  that  at  stage  i there  is  also  a non-empty  set 
of  new  pending  events.  Augmenting  F (R. -SC)  with 
the  part  of  the  resource  increment  flow  network  pertaining 
to  SRi  yields  F (R).  Assume  now  that  F (R.  -SC)  is 
traversed  by  the  flow  resulting  from  flow  contraction.  Even 
if  this  flow  is  maximum  for  F(R.j  - SC),  in  general  it  will 
not  be  maximum  for  F (R)  since  additional  flow  could  be 
pushed  through  edges  with  ee^..  We  call  this  flow 
push-forward  operation  a flow  expansion.  If  at  every  stage 
of  flow  contraction  and  flow  expansion  we  guarantee  flow 
maximality,  we  will  obtain  a maximum  flow  for  F (R)  by 
moving  a minimal  amount  of  flow. 

Flow  Contraction 
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Let  us  call  the  maximum  flow  for  F (RfJ^  In  our 
discussion  we  ignore  the  structure  of  the  flow  sub-network 
for  SCf  by  using  an  auxiliary  flow  network  F^  j that  redirects 
all  flow  entering  dC^  into  the  sink  r.  Formally,  to  obtain  F.  ^ 
we  first  delete  from  F(R)  all  events  in  (5C,.,  together  with  all 
their  incoming  and  outgoing  flow  edges.  We  then  add  an 
auxiliary  edge  e^-^rfor  each  set  of  component  edges 
in  F(RfJ  such  that  e^s  Ri.-SC.  and  e^^SC..  The  capacity  of 
the  auxiliary  edge  e^-^ris  the  sum  of  d\\f^^  . j(ej,  e^)*  We 
call  a function  over  the  edges  of  where ^2) 
is  equal  to/„,^,  if  an  auxiliary  edge,  and 

equal  to  the  edge’s  capacity  if  it  is  an 
auxiliary  edge.  It  is  easy  to  see  that  , is  a maximum 
flow  for  F.  /.  We  call  Res.  ^ the  residual  network  of  F.  j for 
f ' 


For  example,  consider  the  transformation  at  time  t.=4  of 
the  network  in  Figure  4 and  that  at  time  t.^j=3  it 

'^)  = ‘^)  = ej  = I 

arid  Then,  will  not  contain  e^^  and  all 

of  its  incoming  and  outgoing  edged  and  will  contain  two 
additional  flow  edges  of  capacity  c(e^^,  x)  - 3 and  cie^^  x) 
= 7.  The  maximum  flow  function  ^ will  have  the  same 
value  as  for  all  edges  that  already  existed  in  F(R.  ) 
and  for  the  auxiliary  edges  it  is/„,^;  /e^^,  x)  = c{e^^,  x)  - 3 
X)  = c(e^^,  X)  - L Note  that  we  have  only  two 
edges  in  the  residual  network  that  are  associated  to 
auxiliary  edges,  x e,^  with  residual  capacity  3 and  X e^^ 
with  residual  capacity  7,  since  both  auxiliary  edges  are  flow 
saturated. 

We  define  a flow  shift  network  Shift,  as  follows. 


Flow  shift  network:  Shift,  is  a flow  network  with  the  same 
nodes  as  Res^  j.  Shifty  has  a flow  edge  equal  to  a 

corresponding  one  in  Res^j  if  and  Finally, 

for  each  edge  x—^  in  Res^j  such  that  e^x  is  an  auxiliary 
flow  edge  in  Shift,  has  a corresponding  edge  of 
the  same  capacity. 


Shifti  embodies  the  first  route  through  which  flow  is 
pushed  back.  A maximum  flow/^^^/,,;^  , for  represents 
the  maximum  possible  flow  push-back.  After  this  step,  we 
produce  a new  flow  f = for  Fi  t-  Now  we 

need  to  formally  define  another  auxiliary  flow  network, 
Reducci,  to  characterize  the  second  flow  push-back  route, 
the  one  all  the  way  back  to  a.  Let  us  call  Res( Shifti)  the 
residual  network  of  for  f = jf  We  define  a 

flow  reduction  network  Reducei  as  follows. 


Flow  reduction  network:  Reduce^  is  a flow  network  with 
the  same  nodes  as  Res(Shift)  and  edges  ej—^2  identical  to 
Res(Shift)  if  e^x  and  either  €j0fa,  x}  or  Cj=  x and  the 
edge  x—>€  corresponds  to  an  auxiliary  edge  e in 


* Note  that  if  e;  is  a consumer,  an  edge  ej—^XmW  already 
exist  before  the  introduction  of  a corresponding  auxiliary 
edge.  We  tolerate  the  duplication  of  these  edges  since  Fj,j  is 
only  intended  as  a formal  device  for  the  construction  of 
Shifti, 


Using  Shifti  and  Reducei,  we  define  the  Flow.Contraction 
operator  needed  by  the  incremental  envelope  algorithm. 

Flo w.ContractionfF  (RJ , f^^^p  SC^  Aprec ) : 

7 ) Compute  a maximum  flow , Shifty; 

2 ) Compute  a maximum  flow , for  Reduce^ 

3)  Return  f nax.skift.ftfmttx,ndj 

We  now  prove  that  the  operator  keeps  the  flow 
maximum. 

Lemma  4:  The  flow  f maximum  for  £,2. 

Proof:  f is  a flow  of  It  is  also  maximum  since  by 
construction  of  Shifts  it  is  <7,  e)  = 0,  Therefore /Yci; 
-fmax.i-j( o;  e)  and  therefore/’  is  also  maximum  for  £./.□ 

Lemma  is  a flow  for  F(R^^-5C). 

Proof:  is  a flow  for  F._j.  For  it  to  be  a flow  for 

F (Ri  j-Sc)  it  must  be/^„,^/^,  rj=0  if  e—>x  is^  an  auxiliary 
edge.  If  it  were  xj  > 0 for  an  auxiliary  edge,  by 

using  the  flow  conservation  constraint  we  could  show  that 
there  must  be  a path  from  a to  r,  passing  through  e-^x, 
with  all  edges  having  positive  flow.  Therefore,  there  must 
be  a flow-reducing  path  from  r to  cr  in  the  corresponding 
residual  network.  Such  path  is  an  augmenting  path  in  the 
residual  network  of  Reduce^  for  flow  which 

contradicts  the  maximality  of/^^^^^  ..D 

Theorem  6:  is  a maximum  flow  for  F(R.  , - SC. ). 

Proof:  This  is  clearly  true  if/„^^^^,  is  a null  flow  since/’  is 
maximum.  If  is  not  null,  assume  that  is  not 
maximum.  This  yields  an  augmenting  path  77  from  aio  x 
in  F(R,_,  - SC)  forX^„,.  Since jf„.,  is  maximum,  FI  could 
only  have  appeared  after  the  computation  of  f„^,hi/,,r  Since /’ 
is  maximum  for  there  must  be  at  least  one  edge 
belonging  to  FI  that  does  not  belong  to  the  residual  network 
of  for  f otherwise  f„^,hi/,.i  would  not  be  maximum. 
Among  these  edges  consider  the  one  with  either  e^-  x or 
such  that  the  suffix  of  FI  going  from  ^2  to  x has  positive 
residual  capacity  in  Shifti  for  f„^,,,ij,r  A positive  residual  for 
iniplies  that  flow  reduction  pushed  flow  in  the 
opposite  direction,  i.e., ^1)  ^ back-tracing 

e,)  we  find  a positive  flow  path  for  in 

Reducei  from  cr  to  e^.  This  can  only  happen  if  the  capacity 
of  the  path  in  Reducei  is  positive,  which  is  equivalent  to  a 
prefix  path  with  positive  residual  capacity  in  Shifti  for 
fmax,shi/u‘  Tying  the  prefix  and  postfix  at  ^3  yields  an 
augmenting  path  in  Shift,  impossible  since 

is  maximum.^ 

Flow  Expansion 

To  complete  stage  i we  must  now  incorporate  the  event  set 
SRi  to  yield  7?.  and  allow  the  computation  of  P^i  = 

Again,  we  define  an  incremental  operation  on  an 
incremental  residual  flow  network,  the  flow  expansion 
network.  The  network  is  built  on  the  residual  network  of 
F (Ri  i - SC)  for  flow  We  call  this  residual  network 
Res(Contr). 
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Flow  expansion  network:  Expandi  is  a flow  network  with 
nodes  corresponding  to  those  of  Expand^  contains  all 
flow  edges  ej—>e2  in  Res(Contr),  all  flow  edges  in  F (5R) 
and  an  infinite  capacity  edge  e^—^e^  for  each  anti- 
precedence edge  between  Cj  edR^  and  e^  eR^^j-SC^ 

Note  that  by  construction  Expand^  is  the  residual  network 
in  F(R)  for  We  now  define  the  final  operator  needed 
by  the  incremental  envelope  algorithm,  Flow_Expansion. 

Flow_Expansion(F(R^-SC),l^^,  ^,Aprec): 

1)  Compute  a maximum  flow  f^^^^  for  Expand.; 

2) Return 

Theorem  7;/^  , is  maximum  for  F (R). 

Proof:  f^/is  clearly  a flow  for  F(R).  Moreover, z is 
maximum  for  Expandi  and  therefore  there  is  no  augmenting 
path  in  the  corresponding  residual  network.  The  maximality 
of  fmax.i  follows  ffom  the  identity  between  the  residual 
network  of  Expandi  for  f max  exp  i and  the  residual  network  of 
F(Ri)  for  f^^,U 

Flow  Separation  for 

We  can  achieve  further  performance  improvements  by 
minimizing  the  number  of  nodes  and  flow  edges  that  need 
to  be  considered  at  each  stage.  During  stage  i,  two  are 
computed:  after  Fiow^Contractioiij  and  P^.  after 

Flow_Expan$iohj.  We  know  that  each  P^  is  a predecessor 
set  (i.e.,  it  contains  all  of  its  successors  in  the  anti- 
precedence graph),  it  is  flow  isolated  (i.e.,  for  each  pair  of 
events  e,€P^  and  e,)  = 0 mdf^e,.  e,)  = 

0)  and  has  all  exit  edges  saturated  (i.e.,/^(e,  T)  = c(e,  X) 
for  all  eeP^  (Muscettola,  2002).  This  allows  us  to  prove 
that  F can  be  ignored  during  the  computation  of 

FIOW_Contraction^  and  F(P^^^^,J  can  be  ignored  during  the 
computation  of  Flow.Expansion^. 

Let  us  consider  each  maximum  flow  operation  executed 
at  stage  i.  The  first  is  flow  shifting.  Let  us  consider  the 
properties  of  the  set  of  F.  ^ that  contains  the  events  in 
■“  ^ prcdecessor  set  since  all  events  in 

are  at  the  bottom  of  the  anti-precedence  graph  for 
F(R.j).  Moreover,  due  to  added  auxiliary  edges  in  ^ the 
residual  of  the  producers  of  P^  .^j  is  the  same  as  that  in 
and  therefore  equal  to  A(PJrJ).  is  still  flow 
insulated  and  has  all  exit  edges  saturated.  Assume  that  at 
some  point  during  the  flow  shifting  operation  some 
additional  flow  reached  an  event  e*eP^._j.  In  order  for  at 
least  part  of  such  flow  to  reach  x there  must  be  a postfix 
augmenting  path  that  reaches  x from  e\  But  this  is 
impossible  since,  being  P^  ._j  a predecessor  set,  all  postfix 
paths  must  remain  inside  and  all  exit  edges  from 
Emax.i.1  2ire  saturated.  Therefore,  any  maximum  flow 
algorithm  that  searches  for  augmenting  paths  can  avoid 
doing  so  in  Moreover,  in  order  for  any  excess  flow 
pumped  into  events  of  P^  . to  achieve  % that  flow  will  have 
to  be  pushed  back  from  to  Therefore  we  can 

ignore  during  flow  shifting. 


Since  after  flow  shifting  no  flow  has  been  changed  for 
edges  that  touch  P^j,p  P^j.j  maintains  the  flow  insulation 
and  saturation  properties  it  had  before  flow  shifting. 

Considering  now  flow  reduction,  f^„jj  this  can  be 
computed  by  simply  back-tracing  flow  in  F.  j.  Because  of 
the  flow  insulation  of  P^j,p  this  back-tracing  is  either 
performed  exclusively  over  edges  connecting  events  in 
u = ~ o**  is  confined  within  edges 

connecting  events  in  Since  the  entire  flow  that  exits 
P^.,  will  be  back-traced,  the  entire  P^.,  = P^fR.,)  - SC. 

—mox.t-i  ' —max, I-/  max'  i-/'  i 

must  belong  to  the  maximum  predecessor  set 

obtained  after  Flow^Contractioilj.  Therefore  the  contribution 
of  Pmax,i-/  to  can  be  known  in  advance  without 

having  to  modify  any  flow  of  F (P^iJ  during 
Flow^Contraclion..  Hence,  P^(RiJ  can  be  taken  out  of 
consideration  for  future  flow  calculation  as  soon  as  it  is 
computed. 

Finally,  we  can  use  a similar  argument  to  the  one  used 
for  flow  shifting  to  show  that  Flow_Expansion^  can  be 
performed  entirely  over  therefore  allowing  us 

to  ignore  at  any  future  stage. 

Incremental  Computation  of 

We  are  now  ready  to  derive  a recursive  equations  for  the 
incremental  calculation  of  L^naJit)  by  transforming 
Equation  1 through  the  application  of  flow  reduction  and 
expansion. 

Theorem  8:  L^oxfi)  satisfies  this  recursive  equation: 

L^(t)=  A(CJ  + A(P^(R,)) 
iff  = t.  and  i >1 

^max(0  = ^ ^ 

iu 

iv 

ift^b  then 

Emax(l)  = Emax(l”f)» 

Proof:  L^(t)  only  changes  when  R^  changes,  i.e.,  at  a time 
t..  Consider  in  turn  the  application  of  Flow.Contraction^  and 
Flow_Expansion^.  Because  of  flow  separation  after 
Flow_Contraction,„  we  have  = (PJRJ  - dC)  u 

PJ^P  mdPJ  Analogously,  after  Flow_Expansion, 
we  have  P . — P , . uP  (P^  u SR.), 
a)  Fiow^Contraction^:  the  level  after  flow  contraction, 
Fmax.conM  is  the  Weight  of  the  closed  events  after 
confraction  and  of  Pmax.contrj^  Since  Ci  and  Pmax, control  are 

disjoint,  Efnax,contt(li)  ~ ^Ci.j  ty  SCi  (y  P max,contr,i)  “ 

A(Ci,j)  + A(SCi  u(PJiRJ-SC)))  uA(P^JP^JiRJ 
~ SC.)).  Since  for  any  two  sets  A and  B it  is  i4  u(B  - 
A)  - B u(A  - P),  with  B and  (A-B)  being  disjoint  sets, 
we  have  dC^  u (PJRJ  - SCJ  = PJ[RJ 
^SC-PJiK,))-  Hence,  A(Q.,)  + A(dCi  u(PJRJ  - 
SC))  = + A(8CrPJ.RJ).  Since  dC,  cR„  = 

PrJRJ^P'^.JRJ^  it  is  easy  to  see  that  dC,  -P^(R,J 
— SCj  /O  P y^R],).  This  yields  Lmax,conir{^i)  ~ 

^sc,  n I^JRJ)  + A(PJH^J:RJ  - d(C))),  i.e„ 
lines  i,  ii  and  Hi  in  the  theorem’s  statement. 


1CAPS2004  265 


b)  Flow^ExpansioHj:  the  only  new  increment  comes 
from  set  u SR,)  = PJiP^JR,,  - SC,)  u 

SR)  which  yields  line  iv  in  the  theorem’s  statement.^ 


tncremental_Resource„Envelope  (N,  Aprec(N)) 

{ 1:E:={  Group  events  in  the  input  set  N into  entries  Et  with  three 
members:  a time  t and  two  iists  earliest  and  latest  Event 
eeN  is  included  in  Et.earliest  if  et(e)  = t and  in  Eulatest  if 
lt(e)  - 1 Sort  the  Et  in  increasing  order  oft } 

2:  Lmu 0>)  r Maximum  resource  envelope.  V 

3:  tcur  ;=  0;  /*  Current  time  V 

4:  Ud  :=  0;  /*  Envelope  level  at  previous  iteration.  */ 

5:  Lnew  :=  0;  T Envelope  level  at  current  iteration.  V 

6:  Pmax  :=  0;  /*  Maximum  increment  predecessors.  V 

7:Fcur:=0;  /*Resource  increment  flow  graph  with  associated 
maximum  flow  V 

8:  Ecur  :=  0;  r Entry  from  E at  tear.  V 

9:  while  (E  is  not  empty) 

10:  {Ecur  :=  pop(E); 

1 1 : tcur  :=  Ecur.t; 

12:  Lnew  :=  Loid  + weight  (intersection  (Events(F  J,  Ecur.latest)); 

13:  Four  :=  Flow^Conlraction  (Feur,  Ecur.latesL  Aprec(N)); 

14:  <Pmax,  Fcur>  :=  ExUact_P_Max  (Feur); 

15:  Lnew  :=  Lnew  + Weight  (Pmax)r 

16:  Four  :=  Flow.Expansion  (Feur,  Ecur.earliest,  Aprec(N)); 

17:  <Pmax,  Fcur>  :=  Exlract__P_max  (F^^); 

18:  Lnew  :~  Lnew  + Weight  (Pmax)r 
19:  Lmax  I-  append  (Lmaxi  <tcur,  Lnew>)r 

20:  Loid  :=  Lnewp 

) 

return  Lmax? 

} 


Figure  5:  Incremental  envelope  algorithm 


Figure  5 shows  the  pseudocode  of  the  algorithm.  The 
functions  Flow_Contraction  and  Flow^Expansion  receive  as 
arguments  the  current  flow  network  F^^,  which  includes  the 
current  maximum  flow,  the  incremental  set  of  events  that 
need  to  be  added/deleted  Ecur.{earliest,latest},  and  the  anti^ 
precedence  graph  Aprec(N)  for  the  set  of  all  events  N in  the 
plan.  Aprec  carries  the  topological  information  needed  to 
expand  the  flow  network. 

Given  the  current  flow  network  and  its  maximum  flow 
both  stored  in  Four,  ExlracLP^max  returns  both  its 
maximum  increment  predecessor  set  Pmax  and  the  restricted 
network  and  flow  resulting  from  the  elimination  of  the  Pmax. 
Comparing  with  the  formula  for  L^(t)  described  by 
Theorem  8,  line  12  in  the  algorithm  computes  line  15 
adds  Hi  and  line  18  adds  iv.  Note  that  the  pseudocode 
represents  a methodology,  i.e.,  a class  of  algorithms. 
Specific  algorithms  can  be  implemented  selecting  different 
maximum  flow  algorithms  in  Flow_Contraction  and 
Flow^Expansion.  As  we  shall  see  the  worst-case  time 
complexity  of  a specific  instantiation  of  the  methodology  is 
the  same  as  that  the  maximum  flow  algorithm  used. 


Complexity  Analysis 

The  following  complexity  analysis  applies  to  a large 
number  of  maximum  flow  algorithms  used  for 
Flow^Contraction  and  Flow_Expan$ion.  Each  algorithm  has 
a complexity  key,  i.e.,  a measurable  entity  whose  static 
properties  or  dynamic  behavior  during  the  algorithm’s 
computation  determines  the  algorithm’s  time  complexity. 
Table  1 (adapted  from  (Ahuja,  Magnanti  and  Orlin,  1992)) 
reports  the  time  complexity  and  complexity  key  of  several 
maximum  flow  algorithms. 


Algorithm 

Time 

Complexity 

Complexity 

Key 

Labeling 

0{nmU) 

Total 

pushable  flow 

Capacity  scaling 

0(nm 

IorU) 

Total 

pushable  flow 

Successive  shortest  paths 

O(rim) 

Shortest 
distance  to  r 

Generic 

Preflow-push 

0(nm) 

Distance 

label 

FIFO 

reflow-push 

O(n^) 

Distance 
i label 

Table  1:  Complexity  of  maximum  flow  algorithms 


The  Labeling  and  Capacity  Scaling  algorithms  are  based 
on  the  original  Ford-Fulkerson  method.  Their  complexity 
depends  on  the  strict  monotonicity  of  the  flow  pushed 
during  each  algorithm’s  iteration  and  on  the  fact  that  the 
total  pushable  flow  is  bound  by  nil  where  U is  the 
maximum  capacity  of  an  edge  (T—^  or  e—>T.  The 
successive  shortest  paths  class  of  algorithms  is  based  on  the 
original  Edmonds-Karp  algorithm.  The  complexity  depends 
on  the  fact  that  flow  is  pushed  through  augmenting  paths  of 
monotonically  increasing  length.  The  complexity  key  for 
this  class  of  algorithms  is  the  shortest  distance  to  rfor  each 
event  e.  For  these  algorithms  it  is  crucial  to  demonstrate 
that  the  distance  function  d(e)  increases  by  at  least  one  unit 
after  each  iteration. 

Finally,  preflow-push  algorithms  such  as  generic 
preflow-push  and  FIFO  preflow-push  (Goldberg  and 
Tarjan,  1988)  maintain  a distance  labeling  d(e).  These 
algorithms  use  purely  local  operations  that  push  excess 
flow  available  at  node  e^  through  active  edges  e,—>e^  such 
that  d(e)  = diej  + L When  excess  flow  exists  at  some 
node  and  no  edge  is  active,  the  node’s  distance  labeling  is 
increased  by  the  minimum  amount  that  activates  some 
edge.  This  allows  more  flow  to  be  pushed.  The  complexity 
of  preflow-push  depends  on  creating  a valid  labeling  at 
each  iteration  and  on  the  fact  that  for  each  node  the 
distance  labeling  is  monotonically  increasing  up  to  2n-l. 

To  derive  the  complexity  of  our  incremental  envelope 
construction  methodology,  we  analyze  how  each  of  the 
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complexity  keys  described  before  are  modified  across  each 
contraction  and  expansion  stage.  Our  goal  is  to  show  that 
when  moving  across  stages  the  invariant  properties  of  the 
complexity  keys  vary  consistently  to  what  is  required  by 
the  maximum  flow  algorithms  within  each  phase. 

Consider  first  the  cumulative  cost  of  computing  all  flows 
over  2n  stages  respectively  for  SLndfmax,exp,r 

First  note  that  at  each  stage  fmax,red,i  can  be  computed  by 
flow  back-tracing  through  a backwards  depth  first  search 
on  the  resource  increment  flow  network.  Since  this  can  cost 
up  to  0(m),  the  total  cost  of  computing  flow  reduction  is 
0(nm)  and  is  therefore  smaller  than  the  cost  of  applying  a 
regular  maximum  flow  algorithm.  Therefore  we  focus  on 
the  cost  for  the  cumulative  fmax.shift,i  fmax,exp.u  respectively 

^ shift  “ ^ fmax.shiji.i  ^ exp  ~ ^ fmax.exp,v 

Lemma  9:  Neither  nor  are  greater  than  n U. 

Proof:  We  develop  the  argument  for  F^^^  since  the  one  for 
F^^p  is  similar.  Consider  the  total  capacity  of  the  edges 
entering  the  auxiliary  flow  network  Shift..  Its  upper  bound 
is  the  total  capacity  of  edges  with  e~€S(C).  After 
iteration  i all  nodes  in  SC.  are  eliminated  from  further 
consideration,  hence  flow  can  go  through  each  cr-^  in 
Shifty  only  during  iteration  i.  Therefore,  the  total  flow  is 
upper  bounded  by  Z.  \dC\  U = n U.D 

Note  that  the  argument  above  does  not  hold  for  F^^  if  we 
do  not  use  flow  shifting  but  the  flow  is  simply  reduced  and 
then  expanded  again  (Kumar,  2003).  In  this  case  the  same 
flow  could  be  pushed  up  to  n times  and  the  cost  of  F^^ 
would  be  0(nU).  This  would  not  improve  on  the  staged 
envelope  algorithm  in  (Muscettola,  2002). 

Now,  let  us  focus  on  the  second  complexity  key,  the 
shortest  distance  function  d(e)  from  e to  T.  We  will  focus 
on  how  d(e)  changes  at  the  beginning  and  the  end  of 
successive  fmax,shift,i  and  fmax,exp.i  Computations.  Let  us  call 

shift and  ashift/e)  the  distances  at  the  beginning  and  at 
the  end  of  flow  shifting  for  iteration  i.  We  define  similarly 

exp.ii^)  and  df gxp,^^)’ 

Lemma  10:  < (f^/e)  and 

Proof:  Between  the  end  of  flow  expansion  at  iteration  i-1 
and  the  start  of  flow  shifting  at  iteration  the  auxiliary 
flow  network  changes  through  the  elimination  of  F (SC). 
The  elimination  of  edges  could  eliminate  existing  paths  and 
therefore  d(e)  in  the  remaining  residual  capacity  network 
can  only  increase.  Since  Shift,  only  adds  edges  the 
distances  to  r in  Shift,  cannot  decrease  and  therefore 
^ shift  Fo*"  Fxpand^nodQ  distances  to  T could 

further  increase  because  flow  reduction  can  only  eliminate 
residual  network  edges  present  in  Shift,  for  Also, 

from  Lemma  3 the  addition  of  F (^)  cannot  reduce 
distances  since  it  cannot  add  any  edge  from  an  event  in 
Shift,  to  one  in  SR..  Therefore,  < d^^JefU 

Note  that  the  argument  in  Lemma  10  does  not  hold  if 
events  are  added  in  arbitrary  order.  In  this  case  the  addition 


of  edges  can  reduce  the  distance  function  of  some  node  e 
between  a shifting  and  an  expansion  phase.  In  the  worst 
case,  this  may  reduce  some  distance  back  to  a one  unit 
distance  for  each  application  of  maximum  flow  and 
therefore  does  not  improve  on  the  staged  algorithm  in 
(Muscettola,  2002). 

Finally,  consider  reusing  distance  labeling  across 
preflow-pushes  for  shifting  and  expansion,  dexp.i 

and  d^^.  are  the  distance  labelings  at  the  beginning  and  end 
of  shifting  and  expansion.  Assume  also  that  the  distance 
label  of  a node  that  has  not  yet  entered  Expand,  or  Shift,  is 
zero. 

Lemma  11:  j can  be  made  equal  to  d^j.ifor  all  nodes 

in  Shifty  Also,  can  be  made  equal  to  d^Hm 
nodes  in  Expand 

Proof:  The  distance  label  of  a node  remains  valid  when 
edges  are  deleted  or  new  added  edges  always  enter  it  from 
new  nodes.  Also,  at  node  e the  value  of  a distance  function 
d(e)  must  be  an  upper  bound  of  that  of  the  corresponding 
labeling  d(e).  From  Lemma  10  we  know  that  the  distance 
function  can  only  increase  from  Expand^^j  to  Shifts  and  from 
Shifts  to  Expand.,  Therefore  d and  d^^.j^.  are  valid 
choices  respectively  for  and 

We  can  now  prove  the  main  complexity  result. 

Theorem  12:  lncrementaLRB50urce_Envelope  has  time 
complexity  0(Maxflow(n,  m,  IJ))  when  one  of  the 
algorithms  in  Table  1 is  used  for  flow  contraction  and  flow 
expansion, 

Proof:  Assume  we  applied  one  of  the  maximum  flow 
algorithms  in  Table  1 to  the  resource  increment  flow 
network  for  the  entire  flexible  plan  (e.g.,  to  compute  the 
maximum  envelope  level  over  the  entire  time  horizon 
(Muscettola,  2002)).  This  full  maximum  flow  calculation  is 
0(Maxflow(n,  m,  U)).  We  use  Lemmas  9,  10  and  11  to 
prove  that  during  envelope  construction  the  cumulative  cost 
of  using  the  same  algorithm  for  flow  shifting  and  flow 
expansion  is  also  0(Maxflow(n,  m,  U)). 

1.  Labeling  and  Capacity  scaling:  Lemma  9 shows  that  the 

worst  case  bound  for  the  total  flow  moved  during 
shifting  and  expansion  is  at  worst  twice  the  full 
maximum  flow.  Also,  during  shifting  and  expansion 
the  cost  of  finding  an  augmenting  path  is  at  most  m, 
the  same  of  finding  an  augmenting  path  when 
computing  the  full  maximum  flow.  Hence  shifting  and 
expansion  cost  at  most  0(Maxflow(n,  m,  U)). 

2.  Successive  shortest  paths:  the  cost  of  performing  a single 

flow  augmentation  during  full  maximum  flow 
calculation  is  an  upper  bound  of  the  cost  of  a flow 
augmentation  during  shifting  and  expansion.  The 
algorithm’s  complexity  also  depends  on  the  monotonic 
increase  of  the  distance  function  up  to  n after  each 
elementary  operation.  Note  that  until  the  deletion  of  a 
SC.  or  a a node’s  distance  is  bound  by  n,  the  same 
as  during  the  computation  of  the  full  maximum  flow. 
Monotonic  increase  is  guaranteed  by  the  algorithm 
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within  each  shifting  and  expansion  phase  and  by 
Lemma  10  across  phases.  Hence,  the  cost  is 
0(Maxflow(n,  m,  U)), 

3,  Preflow-push  methods:  the  complexity  is  found  through 
amortized  analysis  (Goldberg  and  Tarjan,  1988), 
relying  on  an  appropriate  potential  function  0 and  on 
the  determination  of  its  possible  variations  after  the 
applying  a local  operation  (e.g.,  a saturating  or  a non- 
saturating preflow  push).  One  key  observation  is  the 
monotonic  increase  of  each  node’s  distance  label  for 
each  local  operation.  Both  for  the  incremental  and  for 
the  full  flow  this  increase  is  bound  by  2n-l  and 
Lemma  1 1 guarantees  monotonic  distance  label 
increase  across  phases.  Note  that,  unlike  for  the 
computation  of  the  full  maximum  flow,  for  shifting 
and  expansion  0 increases  also  at  the  beginning  of 
each  shifting  phase,  when  nodes  are  activated  by  the 
creation  of  initial  flow  excesses.  However,  a detailed 
amortized  analysis  (omitted  for  space  limitations) 
shows  that  this  increase  does  not  affect  the  order  of 
complexity  of  the  shifting  and  expansion  phases  that 
remains  0{Maxflow(n,  m,  U)). 

The  worst  case  complexity  of  the  other  phases  of 
lncrementaLRe$ource_Envelope  besides  shifting  and 
expansion  are  dominated  by  0(Maxflow(n,  m,  U).  Flow 
reduction  is  cumulatively  0(nm).  The  total  cost  of 
Extract_P_inax  and  of  incrementally  constructing  and 
deleting  the  flow  network  is  2 0(m),  Finally  the  sorting  of 
events  during  initialization  is  0(n  log  n).n 

Optimized  Preflow-Push  Implementation 

For  preflow-push  implementations  of  the  method,  the 
previous  complexity  analysis  indicates  that  we  need  to 
reuse  the  distance  labeling  function  from  the  end  of  a 
maximum  flow  computation  to  the  start  of  the  next. 

A further  optimization  is  possible.  Consider  the 
maximum  flow  calculation  on  Shift..  During  initialization, 
an  excess  flow  is  loaded  on  each  event  e for  each  edge 
a—>e  in  Shift..  We  know  that  only  a fraction  of  this  excess 
flow  may  reach  r.  The  remainder  will  be  pushed  back  out 
of  Shifty  during  flow  shifting  and  then  pushed  again  through 
the  flow  network  during  flow  reduction.  In  other  words, 
this  flow  travels  twice  through  the  network  before  being 
eliminated.  We  can  remove  this  duplication  as  follows. 
Assume  that,  instead  of  deleting  the  edges  of 
when  constructing  Shift. ^ we  delete  the  edges  of  Shift. 
after  having  performed  the  appropriate  initial  excess 
loading  needed  to  perform  flow  shifting.  In  this  case  the 
flow  that  cannot  be  shifted  will  be  pushed  back  to  the 
source  in  i.e.,  the  source  of  F (RJ,  making  the 
additional  0(nm)  cost  of  flow  reduction  unnecessary. 
Another  possible  optimization  consists  of  combining 
preflow-push  through  Shifts  and  Expand^  by  connecting 
before  running  the  shift/reduce  preflow-push.  In  this  case 
the  initializations  of  contraction  and  expansion  are 
combined  and  a single  preflow-push  is  run  during  phase  i. 


These  optimizations  do  not  affect  asymptotic  complexity 
but  may  have  a significant  effect  in  practice. 

Conclusions 

We  presented  a new  class  of  algorithms  that  efficiently 
compute  resource  envelopes  for  flexible  plans.  The 
methodology  is  0(Maxflow(n,  m,  U))  where  n and  m 
measure  the  size  of  the  activity  plan  and  U measures  the 
maximum  resource  consumption  or  production  of  an  event. 

An  empirical  study  over  all  algorithms  of  the  class  can 
shed  light  on  the  actual  practical  advantage  of  the 
incremental  methodology  versus  the  staged  approach  of 
(Muscettola,  2002)  or  the  incremental  method  in  (Kumar, 
2003).  While  we  expect  that  for  large  problem  sizes  the 
0{n)  cost  reduction  will  be  evident,  practical  improvements 
on  smaller  problems  may  require  careful  design  of  efficient 
and  minimal  data  structures.  Also,  performance  differences 
are  likely  to  occur  among  implementations  that  use 
different  maximum  flow  algorithms. 

By  shedding  more  light  on  the  fine-grain  structure  of 
flexible  plans,  incremental  envelope  construction  is  also 
likely  to  improve  the  “benefit”  part  of  the  cost/benefit 
tradeoff  in  the  use  of  envelopes  for  flexible  scheduling. 

Acknowledgements 

Jeremy  Frank,  David  Rijsman  and  Greg  Dorais  provided 
helpful  comments  on  previous  versions  of  this  paper.  This 
work  was  sponsored  by  the  Automated  Reasoning  element 
of  the  NASA  Intelligent  Systems  program. 

References 

R.K.  Ahuja,  T.L.  Magnati,  J.N.  Odin,  1993.  Network  Flows, 
Prentice  Hall,  NJ. 

R.K.  Ahuja,  M.  Kodialam,  A.K.  Mishra,  J.B.  Orlin,  1997. 

Computational  Investigations  of  Maximum  Flow  Algorithms, 
European  Journal  of  OR,  Vol  97(3). 

T.H.  Cormen,  C.E.  Leiserson,  R.L.  Rivest,  1990.  Introduction  to 
Algorithms.  Cambridge,  MA. 

A.V,  Goldberg,  R.E.  Taijan,  1988.  A New  Approach  to  the 
Maximum-Flow  Problem.  JACM,  Vol.  35(4). 

T.K.S.  Kumar,  2003.  Incremental  Computation  of  Resource- 
Envelopes  in  Producer-Consumer  Models,  Procs.  of  CP2003, 
Kinsale,  Ireland. 

P.  Laborie,  2001.  Algorithms  for  Propagating  Resource 
Constraints  in  AI  Planning  and  Scheduling:  Existing 
Approaches  and  New  Results,  Procs.  ECP  2001,  Spain. 

P.  Morris,  N.  Muscettola,  T.  Vidal,  2001.  Dynamic  Control  of 
Plans  with  Temporal  Uncertainty,  in  Procs.  oflJCAI200I, 
Seattle,  WA. 

N.  Muscettola,  2002.  Computing  the  Envelope  for  Stepwise- 
Constant  Resource  Allocations,  Procs.  of  CP2002,  Ithaca, 

NY. 

N.  Pollicella,  S.F.  Smith,  A.  Cesta,  A.  Oddi,  2003.  Steps  toward 
Computing  Flexible  Schedules,  Procs.  of  Online-2003 
Workshop  at  CP  2003,  Kinsale,  Ireland, 
http://www.cs.ucc.ie/-kbl  l/CP20030niine/onIiticProceedings.pdf 


268  ICAPS  2004 


